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The mean-squared displacement (MSD) of a hard sphere 
and of a dumbbell molecule consisting of two fused hard 
spheres immersed in a dense hard-sphere system is calculated 
within the mode-coupling theory for ideal liquid-glass tran- 
sitions. It is proven that the velocity correlator, which is 
the second time derivative of the MSD, is the negative of a 
completely monotone function for times within the structural- 
relaxation regime. The MSD is found to exhibit a large time 
interval for structural relaxation prior to the onset of the a- 
process which cannot be described by the asymptotic formu- 
las for the mode-coupling-theory-bifurcation dynamics. The 
a-process for molecules with a large elongation is shown to 
exhibit an anomalously wide cross-over interval between the 
end of the von-Schweidler decay and the beginning of normal 
diffusion. The diffusivity of the molecule is predicted to vary 
non-monotonically as function of its elongation. 

PACS numbers: 64.70.Pf, 61.25.Em, 61.20.Lc 



I. INTRODUCTION 

The mean-squared displacement 5r 2 (t) (MSD) is a very 
transparent concept for the discussion of liquid dynam- 
ics jlj. For long times t, Sr 2 (t) increases proportional 
to t and to the diffusion constant of the fluid. In an 
ideal solid, on the other hand, the long-time limit of 
6r 2 (t) is a finite number characterizing the square of 
the particle's localization length. Therefore, the long- 
time behavior of Sr 2 (t) depends sensitively on control 
parameters like density or temperature if the system 
is close to a liquid-glass-transition point. The MSD 
is thus particularly well suited to study glass-transition 
precursors. It can be measured by incoherent inelastic- 
neutron-scattering experiments. However, Sr 2 (t) has to 
be extracted as small-wave-number limit of the inter- 
mediate scattering function [Q, and this makes it very 
difficult to produce accurate data for large time inter- 
vals. Dynamic light-scattering spectroscopy and sample- 
preparation techniques for colloidal suspensions made 
great progress in recent years. It was demonstrated that 
very informative data for the MSD of hard-sphere colloids 
near the glass transition can be obtained ^ ; and promis- 
ing results for this system have also been measured by 
direct-imaging techniques |^,^| . Molecular-dynamics sim- 
ulations are well suited to get accurate data for the MSD 
for liquids near the glass transition. This was demon- 
strated for a binary Lennard- Jones mixture (^] , for a liq- 
uid of diatomic molecules |^|, for models for the van-der- 
Waals liquid orthoterphenyl for the SPC/E model 



for water M, for a hard-sphere-colloid model [Q, and 
for a model of silica ]tT[| . In this paper, general features 
and some quantitative results for the evolution of the 
glassy dynamics as exhibited by the MSD will be consid- 
ered within the mode-coupling theory (MCT) for ideal 
liquid-glass transitions. 

The basic version of MCT is built on approximately de- 
rived closed equations of motion for the auto-correlation 
functions of density fluctuations. The essential input in- 
formations are the equilibrium structure factors which 
are anticipated to vary smoothly with the system's con- 
trol parameters. At certain critical values for the latter, 
there occur bifurcations from solutions for an ergodic liq- 
uid to ones for an amorphous non-ergodic solid. Thus, 
MCT deals with a model for ideal liquid-glass transitions. 
The transition implies a novel dynamical scenario. Its 
features have been worked out by asymptotic solution of 
the equations. The leading-order asymptotic formulas es- 
tablish universal results for the MCT transition like scal- 
ing laws, power-law divergencies for the time scales, and 
anomalous exponents for the spectra |l2],[l3| ■ Many tests 
of the relevance of the MCT results for the explanation 
of the dynamics of glass-forming liquids have been per- 
formed, which are reviewed to some extent in Ref. [p.4| . 
Let us only mention here the recent analysis of data mea- 
sured for propylene carbonate Jig] , studies by means of 
the optical Kerr effect Jl6],fL7| , and the analysis of simula- 
tion data for a binary Lennard- Jones liquid |Tg| ] and for 
silica The outcome of these tests qualifies MCT as 
a candidate for a theory of glassy dynamics and it seems 
justified to continue the preceding studies by unfolding 
some of its implications for the MSD. 

The intention of this paper is to identify further MCT 
results for future tests of this theory. The previous work 
on the hard-sphere system (HSS) |20| — |24f| shall be contin- 
ued by analyzing in detail the MSD for a tagged particle. 
The work on the MCT for molecular systems [^5 will be 
extended by evaluating the MSD for the interaction sites 
of a symmetric dumbbell consisting of two fused hard 
spheres as well as for the molecule's center. The paper 
is structured as follows. In Sec. ||, the equations to be 
solved are listed and the concepts to be used to discuss 
the results are described. Section [nj presents the results 
for the MSD and the analysis of its properties. In Sec. |rv| , 
the findings are summarized. 



II. BASIC FORMULAS 
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A. Description of the system 

In this section, the systems to be studied and the func- 
tions to be used for a description of their dynamics shall 
be defined. A system of N atoms of mass m distributed 
with density p is considered as solvent. The points in 
configuration space are specified by the particle posi- 
tions F K , k = 1,2, • • • , N. The basic variables for the 
description of the structure are the density fluctuations 
with wave vector q, pq = X) K ex P(*9' ^k)- ^ (' ") °-C- 
notes canonical averaging for temperature T, the struc- 
ture factor is S q = (| p$ \ 2 )/N, where q =| q | is the 
wave number. As simplest example for a solute, a tagged 
atom of mass ms and position F$ shall be considered. 
The distribution of the atom is described by the den- 
sity fluctuation p~ = exp(iq ■ fg). The solute-solvent 
interaction shall be characterized by the direct correla- 
tion function c q = (/?§* pq) /(pS q ). As a more compli- 
cated solute, a symmetric rigid diatomic molecule shall 
be chosen. Its position is specified by two interaction 
sites (Fa, rs), which have the same mass ttia- The posi- 
tion of the molecule can also be described by its center 
Fc = (Fa + Fb)/2 and the unit vector e= (Fa — Fb)/L, 
where L = \ Fa — Fb \ ■ The configuration variables for 
the molecule can be built with the two density fluctu- 
ations pi — exp(i<f • F a ); a — A, B. Equivalently, one 

N 



can use the number fluctuations pi 



(p4 + p$ )A/2 

and the "charge" fluctuations p| = (/>4 - p?)/y/2. The 
solute-solvent interaction can be characterized by the 
direct interaction-site-solvent correlations c q = c q = 
(P$*Pit)/\pS q v$] @. Here w»> z = l±sm(qL)/qL de- 
note the intra-molecular structure factors. 

Three kinds of mean-squared displacement functions 
of time t shall be discussed, Sr 2 (t) — ([F x (t) — ^(O)] 2 ). 
Here and in the following, the label x = S, C and A refers 
to the position of a tagged particle, of the center of the 
molecule, and of the atomic center in the molecule, re- 
spectively. It will be more convenient to use the following 
abbreviation 



The non-trivial time dependence of F x (t) — F x (0) comes 
about since the forces on the solute fluctuate in time; and 
this is caused by the density fluctuations of the solvent 
and by the fluctuations of the probability density of the 
solute constituents. These quantities are described by 
the density correlator <j> q (t) = (Pq(t)* Pq) / (\ Pq | 2 ) of the 
solvent, by the tagged-particle-density correlator (j> q (t) = 



(pf(t)*Pg), and by the molecule's correlators 



kN,Z 



(*) 



,N,Z 



\Pq 0) Pi[ )/ 

the desired functions, since []l|,t!5| 



These correlators also determine 



b S q (t) = l-q 2 A s (t) + 0(q i ), 



l-q 2 A c (t) 



0(q A 



^(t) = C 1 (t) + 0(q 2 ). 



(5a) 
(5b) 
(5c) 



M*) = l<fr 2 (f), x = S,C, A. 



(1) 



B. MCT approximations 

In this section, those equations shall be listed which 
have to be solved numerically. Within the Zwanzig-Mori 
theory, an exact equation of motion can be formulated for 
the density correlator: df(j)q(t) + Vt 2 q tf> q (t) + J* * dt' M q (t - 
t') d t '4> q (t') = 0. Here VI 2 = v 2 q 2 /S q with v denoting the 
thermal velocity specifies a characteristic frequency f2 g , 
and M q (t) denotes a fluctuating- force correlator called 
the relaxation kernel Q). Within MCT, the kernel is split 
into a regular part M q es (t) dealing with normal-liquid 
effects and a mode-coupling kernel Q 2 m q (t) describing 
the cage effect. If one introduces an operator TZ for the 
regular dynamics by K4> q (t) = [df<p q (t) + J* dt' M™ B (t - 
t') dt<<j)q(t')]/£l 2 , one can write 

+ + / dt> ~ Bt'M 1 ?) = °- ( 6 ) 
Jo 

The crucial step in the derivation is the application of 
Kawasaki's factorization approximations to express the 
kernel m q (t) as mode-coupling functional T q of the cor- 
relators 



The MSD of the molecule's constituents can be decom- 
posed into one contribution due to translation of the cen- 
ter and one due to reorientation of the axis |E5l 



where 



A A (t) = A c (t) + ±L 2 [l-Ci(i)], 



Ci(i) = (e(*)-e) 



(2) 



(3) 



is the dipole correlator. The time derivatives of the MSD 
provide the velocity-auto-correlation function Q . Let us 
consider the one for the velocity vs(t) of the tagged par- 
ticle only, Ks(t) — (Fs(t) ■ Vs), where 



d?A s (t) = ±K s (t). 



(4) 



m q (t) = T q [(j)k(i)\ 



2(2tt) s 



dkV(q;k,p) fkfp 



(7a) 
(7b) 



Here p abbreviates q — k. The coefficients V(q; k,p) are 
given in terms of the structure factor |^7j . 

None of the MCT results for structural relaxations, 
in particular none of the universal results to be cited in 
the following Sec. [IC, depend on the model for M q ° s (t). 



This holds up to the choice of some time scale to be 
denoted below as to. But, the kernel M q ° s (t) shall be 
specified in order to have controllable quantitative results 
for all times. Specifically, a model with M q C6 (t) = shall 
be chosen. The operator 1Z shall be complemented by 
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an index H indicating that a Hamiltonian dynamics is 
considered for the short-time motion 



n n 4> q (t) = d^ q (t)/n 2 q . 



(8a) 



This model overemphasizes oscillation features. A more 
realistic model would include at least some friction term 
as it is caused for low frequency phenomena by binary 
collision events. But, no detailed proposals for a treat- 
ment of such effects have been made so far within MCT. 
Some results also will be presented for a simplified colloid 
model. Here, the inertia term from Eq. ( pa| ) is neglected 
and the regular term is chosen as a ^-independent white 
noise kernel. It is explained in more detail in Appendix ^ 
that this model corresponds to the conventional treat- 
ment of colloids by coarse-graining the time over inter- 
vals of the duration of collisions of the solvent molecules 
with the mesoscopic colloid particles. As a result, the 
short-time motion is treated by a Brownian dynamics 



n B ^> q (t) = Tq d t ct> q {t). 



(8b) 



Here r q = S q /(D q 2 ) with D denoting the single- 
particle-diffusion constant. 

Equations (g) and (g) hold analogously for the so- 
lute correlators (j) q (t), x — S, N and Z. One gets 
ilq 2 = v 2 s q 2 with vs denoting the tagged particle thermal 
velocity. The relaxation time for the Brownian motion 



is r, 



s _ 



1/(Z)q<7 2 ) with Dq denoting the tagged parti- 
cle short-time diffusivity. The more involved expressions 
for the characteristic frequencies Sl q ' Z can be found in 
Ref . |2[| . Brownian dynamics shall not be considered for 
the dumbbell molecule. The fluctuating-force kernels are 
functionals of the correlators 4> q (t) and 4> q (t): 



1 



^nifkJp 



(2tt) s 



dkV x {q;kJ)fUp 



(9a) 
(9b) 



for x ~ S,N, Z. Again, p abbreviates q — k and the coef- 
ficients V x (q; k,p) are given by S q and the direct correla- 
tion functions |2^j2^]. It is cumbersome to calculate the 
required q — > limits in Eqs. (|^) numerically from nu- 
merical solutions for 4> q (t)- It is more adequate to carry 
out the limit analytically in the equations of motion for 
4> q (t) so that one gets equations of motion for the de- 
sired functions. The non-trivial parts of these equations 
are convolution integrals defined with the q — > limits of 
the kernels m x (t) . One gets for the dipole correlator 

d 2 C 1 (t) + 2v 2 R C 1 (t) 



+ 2vl / dt'mz(t-t')d t <Ci(t') =0, (10) 
Jo 

where wr is the thermal angular velocity of the 
molecule Q . From the equation for c/) q (t) one get s a 
Zwanzig-Mori equation for the velocity correlator 



d t K s (t) + v 2 s / dt'm s (t-t')K s {t') = 0. (11) 
Jo 

Integrating twice over t one gets with the aid of Eq. (jij) 
and the initial conditions As(0) = and 9 t As(0) = 0, 

H: d t A s (t)-v 2 t + v 2 f dt'ms(t-t')A s (t') = 0. 

Jo 

(12a) 

The corresponding equation for Brownian short-time dy- 
namics p3| is derived in Appendix ^| 

B: A s (t)-D$t + D$ f dt'm s (t-t')As(t') = 0. 
Jo 

(12b) 

The same procedure leads to the equation of motion for 
the MSD of the center 

H: d t A c {t) -v^t + v 2 , [ dt'm N (t-t') A c {?) = 0, 

Jo 

(13) 

where the kernel is denoted by m.]\[(t) and vt is the 
thermal velocity for the molecule's translation p5j. The 
problem of the molecule's dynamics reduces to that of 
a tagged atom with wis — 2m a, if the limit L — > is 
considered. 

The kernels in the preceding Eqs. (|Io|)-(|l3|) are given 
by mode-coupling functionals 



m x {t)=F x [<j>t{t),j> v {t)\, 



(14a) 



? x [ft , fp\ = ^2 / o dk fc4 P S >< vX ( fc ) fk h > (14b) 

for a; = S,N,Z. Rerev s (k) = cf 2 ^,v N (k) = 2c£ 2 w% 
and v z {k) = (L 2 /6)c£ 2 w z H§. A fluctuating force with 
vanishing wave vector can couple to density fluctuations 
of the solvent for all wave vectors k provided the atom 
or molecule can absorb the recoil with wave vector — k. 
Therefore, one needs the superposition of density cor- 
relators 4>k{t) and (j)%{t) for all wave numbers k for the 
calculation of the kernels m x (t). 



C. Universal results 

Universal properties of the MCT-glass-transition sce- 
nario are formulated by the leading-order asymptotic ex- 
pressions for the long-time dynamics for states near the 
transition points. This paper focuses on features beyond 
the universal ones, but the universal formulas shall be 
used as reference. In this section, those formulas |12 j | sh all 
be compiled which are needed in the following SecT|l| for 
the description of the results. 
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The equilibrium structure of the system may depend 
on, say n, control parameters which can be combined to 
a control-parameter vector V. A separation parameter 
cr(V), a smooth function of V, can be defined with the 
aid of the mode-coupling functional T q . For states with 
control parameters V such that a < 0, correlation func- 
tions decay to zero: 4> q (t — > oo) = 0. But for states with 
<r > 0, density fluctuations exhibit spontaneous arrest: 
<j) q {t — > oo) = f q > 0. The Debye- Waller factor / g is 
to be evaluated from the mode-coupling functional T q in 
Eqs. (0) via the equation f q /(l - f q ) = T q [f k ] ©■ The 
set of critical points V°, defined by cr{V c ) = 0, separates 
liquid states from glass states. This result holds for all 
correlators cf) A (t) = (A(t)* A) / (\A\ 2 ) of variables A cou- 
pling to density fluctuations. While <f»j\.(t — > oo) vanishes 
for the liquid, generically, the limit f A = 4>A(t - > oo) is 
positive for the glass. If A refers to pg, p| or e, $a de- 
notes the Debye- Waller factor f q , the Lamb-Mossbauer 
factor fy or the Edwards- Anderson parameter fx, re- 
spectively, of the glass. Crossing the transition points, 
the long-time limit changes discontinuously from zero to 
the critical value f c A . For the variation of f A upon ap- 
proaching the transition from the glass side, one gets for 
small a in leading order 



fA = fl + h A 



(15) 



Here hA is called the critical amplitude for variable 
A. Every point V° is characterized by a number A, 
< A < 1, which is called the exponent parameter. 
The significance of A is explained below in connection 
with Eqs. (pj)-( jl8l). It is a matter of convention not 
to incorporate \/l — A in the amplitude hA, in order to 
simplify some of the following formulas. The quantities 
f A , hA and A are calculated from T q . They are equilib- 
rium quantities which are the same for Hamiltonian and 
Brownian dynamics. 

The parameter A determines an anomalous exponent 
a, < a < 1/2, which is called the critical exponent. 
There holds the equation T(l — a) 2 /r(l — 2a) = A, where 
r denotes the gamma function. The long-time decay at 
the critical point is given, up to corrections of order t~ 2a , 
by the power law 

Mt) = fl + h A {h/t) a , (7 = 0, t/to » 1. (16) 

The A-independent time to is the relevant microscopic 
scale for the bifurcation dynamics. It depends on all 
details of the transient dynamics as well as on the mode- 
coupling functionals for parameters at the transition 
point. 

The first scaling law of MCT deals with the dynamics 
for small a in a time interval where r\ = 4>A(t) — f A is 



small. In a leading expansion for a 
gets 



and r\ — > 0, one 



Mt)=f A + h A G(t), 
G(t)=y/Wg±(t/t a ), a^0, 
ta=t Q /\a\ s , <5 = l/2o. 



(17a) 
(17b) 
(17c) 



The functions g±(t) are determined by A. Thus, the 
control-parameter dependence of the dynamics is deter- 
mined by the correlation scale yf\tJ\ and by the first crit- 
ical time scale t a . One gets g±(t — * 0) = l/t a , so that 
Eq. ( |l6| ) is reproduced for fixed large t if a tends to zero. 
Since g + (i — > oo) = — A, also Eq. ( p"5| ) is repro- 

duced. 

There holds g_(t -> oo) = -Bt b + 0(1 /i b ). The 
anomalous exponent b, < b < 1, which is called the 
von Schweidler exponent, is to be calculated from the 
equation T(l + b) 2 /T(l + 2b) = A. The constant B is 
of order unity. Substituting this result into Eqs. (|Tt| ) 
one gets von Schweidler's law for the decay of the liquid 
correlator below the plateau f\ 



<P A {t) = f c A -h A (t/t' a ) b , t CT «i, 



-0. (18a) 



The control-parameter dependence is described by the 
second critical time scale 

t' a =t Q B- 1 ' b /\a\\ 7 =(l/2o) + (l/26). (18b) 

Following the terminology of the glass-transition liter- 
ature, the decay of 4>A(t) below the plateau f A is called 
the a-process. For this process, there holds the second 
scaling law of MCT in leading order for a — > —0: 



v (t) = 4> A (i), t = t/t' (T , t a <$:t. 



(19) 



The control-parameter-independent shape function <pA (t) 
is to be evaluated from the mode-coupling functionals at 
the critical points V c . The differences of the dynamics 
as they are caused by different models for the short-time 
dynamics merely enter via differences in the scale to . For 
short rescaled times t, one gets (f) A (t) = fl ~ + 
h' A t 2b + • • •, so that Eq. (18a) is reproduced. The ranges 
of applicability of the first and the second scaling laws 
overlap; both scaling laws imply von Schweidler's law for 

t a < t < t' a . 

Suppose, the system is driven through the transition 
point V c by smooth variation of some parameters 9 like 
the temperature, the density or, for a colloid, the salt 
concentration of the solvent. Let 9 C denote the value 
where V(9 C ) — V c . Then one can write for a, in leading 
order for small (6 — 9 c ), the expression a = Cg (9 — 9 c )/9 c . 
The constant Cg depends on the choice of 9 and connects 
the distance parameter e — (9 ~ 9 C )/9 C with the relevant 
separation parameter a. 



D. The model 

Hard spheres of diameter d shall be used as model for 
the solvent atoms. For this case, all equilibrium quan- 
tities are specified by the packing fraction ip — irpd 3 /6. 
The MCT model for the hard-sphere system (HSS) shall 
be defined by two further technical assumptions. First, 
the structure factor S q and the direct correlation func- 
tion c q are evaluated within the Percus-Yevick theory . 
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Second, the wave numbers are discretized to 100 equally 
spaced values qd = 0.2, 0.6, 1.0, • • • , 39.8. The details of 



the transformation of the functional in Eq. (7b) to a poly- 
nomial in the 100 variables <f> q (t) can be found in Rcf. p2| . 
Representative solutions are shown in Refs. J2J] and |2q | 
for the Hamiltonian dynamics and in Ref. |22j| for the 
Brownian dynamics. There is a liquid-glass transition at 
the critical packing fraction tp c w 0.516 |2^,|27j]. For the 
exponent parameter one gets A = 0.735, and this implies 

a = 0.312, b = 0.583, B = 0.836, S = 1.60, 7 = 2.46. 

(20a) 

For the separation parameter a one gets in leading order 



a = 1.54 e, e = (ip - tp c )/<Pc 



(20b) 



The microscopic time scales t for Hamiltonian dynam- 
ics |2?| and Brownian dynamics 22 1, respectively, are 



t% = 0.0236 (d/v), = 0.00265 d 2 /D . 



(21) 



As atomic solute, a tagged particle of the solvent shall 
be considered, i.e., ms — m and Cg is identical with the 
direct correlation function c q of the HSS. As molecular 
solute, a symmetric dumbbell of two fused hard spheres of 
diameter d and mass tua = Tn shall be chosen. The elon- 
gation Q — L/d shall be used as control parameter of the 
solute. The wave numbers are chosen discrete as above. 
The direct correlation function is expressed as series of 
contributions ce(q), i — 0, 1, • • •, obtained by expanding 
the molecule-solvent correlations in spherical harmon- 
ics p5| . The sum over I is truncated at £ co = 8. The 
ce(q) are evaluated with the Percus-Yevick theory P9[ . 
Representative results for molecule's correlators cf>g (t), 
tfiq (t) and C\ (t) are shown in Ref_j£5[ , and for </>^ (t) for 
the Brownian dynamics in Ref. |23[ . 

The figures to be discussed below and the numbers to 
be mentioned are evaluated for the above-specified model 
as follows. First, for a representative set of packing frac- 
tions, Eqs. (||) and (Q) are solved for the density cor- 
relators <p q (t), and this for Hamiltonian as well as for 
Brownian dynamics. These correlators are used to de- 
fine the kernels in Eqs. ([)]), so that, as second step, the 
tagged particle correlators 4>g (t) could be evaluated, also 
for both examples for the short-time dynamics. Further- 
more, for every tp, the equations for molecule's correla- 
tors 4>q' z {t) are solved for 10 values for the elongation 
(. These results are substituted into Eqs. ( |l4|) for the 
kernels m x (t) so that, as last step, Eqs. ([l0|)-([13|) for the 
desired functions A x (t), C\(t) and Ks(t) can be solved. 



III. RESULTS 

A. The diffusion-localization transition 

If a tagged particle would experience a mere Newto- 
nian friction force, the velocity correlations would decay 



exponentially, Ks(t) cx exp[— (t/r)]. The cage effect in 
dense liquids manifests itself by a qualitatively different 
behavior, namely by oscillatory variations with a decay 
of Ks{t) to negative values M. Figure |l](a) demonstrates 
this phenomenon. With increasing density, the crossover 
time to negative values shortens and the damping of 
the oscillations increases. A Green-Kubo formula relates 
the particle diffusivity D$ to the zero-frequency velocity 
spectrum: D$ — (1/3) J °° dt Kgit) Q]. Negative contri- 
butions to Kg(t) reduce the diffusivity with increasing 
tp. From Eq. (pi]), one gets Dg as inverse of the zero- 
frequency spectrum of the relaxation kernel fll2|| : Ds = 
l/f£°dtms(t). From Eqs. @, one obtains for the long- 
time asymptote of the MSD: limt_ >0 o &s(t)/t = Ds 0]. 

For glass states, density fluctuations arrest for long 
times: 4>g(t — > 00) — fg > 0. The Lamb-Mossbauer fac- 
tor fjj 1 is to be evaluated from the mode-coupling func- 
tional J-q in Eqs. (g) via the equation fg /(l — fg) — 

^qifk'fp] @- It approaches unity for q tending to 
zero. A localization length rg can be introduced to 

fS_ 



characterize the width of the -versus-?? curve: fg 

{qr s ) 2 + 0(q 4 ). From Eq. Q, one gets for the MSD: 
2 



t — >oo 



■ s- 



Using Eqs. (|12|), one can express 



1 

lim 

r| as inverse of the long-time limit of the relaxation ker- 
nel @: r| = l/m s (t -»• 00) = l/^ s [/f , /,]. If the den- 
sity increases, the localization length r$ decreases. As a 
result, the frequency of the oscillations of the particles 
in their frozen cages increases, as is shown in Fig. ^(b). 
But, in contrast to what is found for liquid states, the 
damping of the oscillations decreases upon compression. 
This reflects the formation of anomalous oscillation peaks 
in the density-fluctuation spectra, which have properties 
of the so-called boson peaks of liquids and glasses Q . 

The ideal liquid-glass transition implies a transition 
from a regime with particle diffusion for ip < (p c to one 
with particle localization for tp > tp c . The former is 
characterized by Ds > and 1/rg = and the latter 
by Ds = and 1/rs > 0. The subtleties of the glass- 
transition dynamics occur outside the transient regime. 
They can be discussed best on logarithmic scales as in 
Fig. ||. For very short times, say t < to, interaction 
effects are unimportant and limt_>o As(t)/t 2 — vf/2 re- 
flects ballistic motion. For times larger than to, the cage 
effect leads to a suppression of As(t) below the short- 
time asymptote. For such large times that 5r s (t)/d 2 
reaches unity, the MSD approaches the diffusion asymp- 
tote, A,g(i) w Dst, as is shown by the dotted straight 
lines drawn for the curves with labels n — 1 and n = 9. 
Upon increasing ip towards tp c , the diffusivity decreases 
towards zero. Figure || shows that there holds the power 
law D s h cx I e I for I e I < 0.1. 

The lowest line in Fig. || deals with the same glass 
state <p = l.lp c , which was considered in Fig. |l](b) for 
the label n — 3. For this density, there is no obvious 
glassy dynamics. Rather, As(t) has approached its long- 
time limit r| after the oscillations have disappeared for 
t rs 1. Decreasing ip towards <p c , the softening of the glass 
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manifests itself by an increase of the localization length 
rs- At the transition point tp = p c , the critical value 
r c s = 0.0746 d is reached. This upper limit for rs is con- 
sistent with Lindemann's melting criterion j27j. Using 
Eq. (|l5|) for f q and f q and substituting these formulas 

into F s [fk i fp]> it follows that the glass instability at ip c 
causes a •y/a'-anomaly for the localization length, 



l-A + 0(tr), 



(22) 



where r c s 2 = 0.00557 d 2 and h s = 0.0116 d 2 . Figure | 
demonstrates that the leading asymptotic formula ac- 
counts for the Tg-veisns-ip dependence for e < 0.01. But 
the data for n < 4, i.e., for e > 0.05, are not anymore 
described by the y/a— law. The range of applicability for 
the asymptotic description of r 2 s is remarkably smaller 
than the one for the corresponding description of f q for 
intermediate wave numbers p2] | . 

The glass curve for e = 0.01, shown in Fig. ||with label 
n = 6, exhibits a decay between the end of the transient 
oscillations and the arrest at r 2 s which is stretched over a 
time interval of about two orders of magnitude. A similar 
two-decade interval is needed for the liquid curve with 
label n = 6 to reach the critical value r| 2 . After crossing 
r s 2 , two further decades of an upward bent log 10 Ag(i)- 
versus-log 10 t variation is exhibited before the diffusion 
asymptote is reached. The indicated slow and stretched 
time variation is referred to as glassy dynamics. 



B. The structural-relaxation regime 

For times outside the transient regime, say t > C*to, 
the density correlators can be written in the form 4> q (i) — 
^o(*/*o)- Here t is the scale introduced in Eq. (16). 



The functions <fi* are determined uniquely by the mode- 
coupling functional J- q , i.e., they are given by the equi- 
librium structure. This holds for all choices of the reg- 
ular kernels in Eq. (^|), in particular for the two models 
defined by Eqs. (g) and @ ||||2|]3(]fl|] ■ Correspond- 
ing results hold for the density correlators of the tagged 
atom |30) and of the molecule |^5| . The solutions of the 
specified MCT model for colloids arc completely mono- 
tone p2j. A function F(t), defined for t > 0, is called 
completely monotone if all derivatives exist and 



{-d/dt) n F{t) > 0, n = 0, 1, 



(23a) 



According to Bernstein's theorem this property is 
equivalent to the existence of a distribution ,0(7) > so 
that 



F(t) 



e 7t p(j) d*y. 



(23b) 



Thus, one can write for t > C*to: 4>q{t) — 
Io° e -7 ^* - 1 Pqij) d'y with p q (j) > 0, i.e., the functions 
(f>* deal with relaxation. Corresponding formulas hold for 



the solute correlators cj) q (t) . Representative examples for 
Pqij) are discussed in Ref. Q. 

From Eqs. (§), one gets d{t) = C{{t/t Q ) for t>C*t , 
where C* is completely monotone, and also 



A x (t) = A* x (t/t ), t>C*t Q , 



S, C, A. 



(24) 



Here, the functions C{ and A* are determined by the 
equilibrium structure. With Eq. (g) one can express 
the velocity correlator in terms of the structure function 
F(t) = -d 2 A s (r): 



K s (t) = -3F(t/t )/t 2 , t>C*t . 



(25) 



It can be shown that the function F is completely mono- 
tone. The proof does not provide further insight and is 
delegated to Appendix |X| 

In Ref. p8| , the curves for <p q (t) have been compared 
with the ones for cj) q (t/to). The interval for structural 
relaxation was larger for the colloid model than that for 
the model with the underlying Hamiltonian dynamics. 
Thus, the structural relaxation interval can be identified 
as the one, where the specified curves for the two models 
collapse. Figure g exhibits such comparison for the MSD. 
For I e I < 0.01, the curves agree within the accuracy of 
the drawing. This holds, provided t > 20t , i.e., C* ~ 20. 
The result is nearly valid also for distance parameters as 
large as | e | = 0.1. But there is a small offset between 
the full and the dashed liquid curves for x = 1. This 
means, that there is a smooth drift of to with changes of 
p, which is different for the hard-sphere colloid and for 
the conventional HSS. 

In Fig. H the rescaled velocity correlator £q Ks(t)/3 for 
the critical packing fraction is shown as full line. This 
diagram is an extension and magnification of the dotted 
lines from Fig. for t > 0.2. The dashed line is the 
analogous result tg 2 d 2 A c s ol (t) calculated for the colloid 
model. The latter function is (t^D§) 2 F col (t), where the 
completely monotone function F col (t) was in troduc ed i n 
Appendix [X] in connection with Eqs. (A8a) and flAS|). 
According to Eq. ( p5|) , the two curves should collapse on 
the function F(t/to). This is the case for t > 20to within 
small error margins. The curves demonstrate stretched 
relaxation to zero, which cannot be adequately repre- 
sented on linear scales. It is shown also that oscilla- 
tory motion tends to mask glassy relaxation. From now 
on, the discussion will focus on the structural-relaxation 
regime t > 20to- 



C. Scaling- law descriptions 



In this section, 



it shall be examined how well the 

can ac- 



lcading-order asymptotic results from Sec. II C 



count quantitatively for the MSD. Let us start with As(i) 
for <p — p c . This result for the critical dynamics, i.e. the 
dotted line in Fig. ^, is reproduced as full line in the 
semi-logarithmic presentation in Fig. ^. The transient 
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dynamics for t < 20io accounts for about 45% of the total 
increase of Ag(i) from zero to the long-time asymptote 



„c 2 



c 2 



The structural relaxation needed to approach r| 
up to 5%, i.e. 50% of the total increase, is stretched over 
a large interval of about four orders of magnitude time 
variation. The leading-order formula for the MSD at the 
transition point is in analogy to Eq. (|l6|): 



r S 



h s {t /t) a , tp = <p c , t»*o- (26) 



The dashed line demonstrates that Eq. ( |2q ) describes 
- within a 5% error margin - about 25% of the to- 
tal increase of Ag(t). There remains the large part of 
the structural-relaxation interval between 20to and t* — 
792io which is not adequately accounted for. About half 
of the structural-relaxation increase of Ag(i) is outside 
the range of applicability of the leading-order asymptotic 
formula. Extending Eq. ( p6| ) by inclusion of the leading 
correction, one gets a description of the MSD up to errors 
of order (W0 3q : A s (t) = r% 2 - h s (t /t) a + k s (t a /t) 2a . 
From Ref. [§3| one deduces k s = 0.0143 d 2 . The dashed- 
dotted line shows how inclusion of the correction term 
expands the range of the analytic description. 

The first scaling law for the MCT-bifurcation dynam- 
ics, Eqs. (|l7j), implies with Eq. ( [5a| ) 



As(t) 



a < 1 



a \g±(t/t<r), 
t > t . 



(27) 



For the liquid states with a < 0, it describes how the 
As(£)-versus-£ curve crosses the plateau r c s 2 . For the 
glass states with a > 0, it describes the approach towards 
the arrest at r| = As(t — > oo). The control-parameter 
independent functions g±(t) for the HSS value of A are 
discussed in Fig. 10 of Ref. |2y] . The dashed lines in Fig. |?j 
exhibit Eq. ( p7[ ) for three liquid states. They agree with 
the MCT solutions within a 5% error margin within the 
intervals marked by diamonds. Within these intervals, 
A s (t) increases from about 0.0043 d 2 to about 0.0086 d 2 . 
Formulas like Eqs. ( |l6| ) and (|26| ) are the basis for the 
derivation of the MCT-scaling laws. Therefore, it follows 
from Fig. ^ that the part of the structural relaxation 
regime between 20to and t* remains outside the range of 
validity of Eq. @>. 

For small rescaled times t = t/t a , Eq. (|27| ) reproduces 
Eq. ( p6| ) for o — > 0. The r.h.s. of Eq. Q27|) becomes in- 
dependent of a and agrees with Eq. (|2^). This explains, 
why the dashed lines in Fig. ^ for x = 3 and 4 collapse for 
t < 10 and why the corresponding diamonds are located 
near t* . For x = 2, the separation parameter a is already 
so large, that g±(t* /t a ) differs remarkably from (t rT /t*) a . 
Therefore, As(t) does not reach the t~ a asymptote for 
t w i* , and the corresponding diamond shifts away from 
t*. The x = 3 curve follows the critical asymptote from 
Eq. (^) for a time interval of less than two decades. Such 
interval would not be large enough for a compelling ex- 
perimental confirmation of the t~ a law. To identify the 
t~ a law in its pure form for the model under study, | e | 
must not exceed 10~ 4 . 



For large rescaled times t, the master function for the 
glass approaches 9+(i — > oo) = 1/Vl — A, and Eq. (|27| ) 
reproduces Eq. (p2). According to the preceding para- 
graphs, this explains the solutions for the glass provided 
the long-time limit r| is located in the shaded interval of 
Fig. 0; and this is demonstrated in Fig. |[ 

Function g-(t) for the liquid is zero for t- — 0.704. 
Thus, As(t-t a ) = r c s ? and the interval for the increase 
of Ag(t) to the plateau value r^ 2 expands proportional 
to t a if 1^9 increases towards (f c . For large t, one gets 
von-Schweidler's law 



As(t) 



r c2 
r 5 



(t/O 



-0, to- < t < t' 



(28) 



Therefore, the long-time end of the range of applicability 
of Eq. pTl) expands proportional to t' , as is indicated 



by the filled squares in Fig. 0. Formula ( p8| ) is exhibited 
by the crosses. These approach the plateau r c s 2 for t <C 
t CT and the dashed scaling-law lines for t S> t a . Since 
t' a /t a — > oo for | e | — ► 0, the time interval for the von- 
Schweidler-la w de scription expands with decreasing | e |. 

Equation (|l7a| ) formulates the factorization theorem 
for 8cj>A(t) = 4>A{t) ~ Ia : m a leading-order expansion for 
small S<p a, the deviation 54>A(t) of the correlator from 
the plateau value j c A factorizes in a control-parameter- 
independent amplitude Ha and a function G(t). The 
function G(t) is the same for all variables A and describes 
the time-and-contr ol-pa rameter dependence of <f>A(t) by 
a scaling law, Eq. (|l7b|) . This theorem can be tested by 
identifying the time interval and the range of distance pa- 
rameters e for which the diagrams for <f>A{t) = S<t>A(t)/fiA 
collapse with G(t). Figure H demonstrates such a test for 
the dipole correlators C\ it) for three values of the elon- 
gation parameter £. The markers for 20io, t* , t a and 
t' a have been added to ease a comparison with Fig. [t]. 
Obviously, the scenario for the plateau crossing is the 
same for C\(t) as discussed above for As(t). To corrob- 
orate this conclusion, the rescaled result for the MSD, 
A s (t) = [r c s 2 -A s (t)]/h s fore = -0.001, has been added 
to the figure. The x = 2 results show that for negative 
Ci(t) the full lines follow the sequence £ = 0.6, 0.8 and 
1.0 from top to bottom, where the latter two curves are 
very close to each other. The same behavior is observed 
for positive C\(i). This observation exemplifies a general 
implication of the leading corrections to the factorization 
theorem 

A possibility for the definition of a characteristic time 
scale ta for the a-process of variable A is given by the 
time needed to complete 95% of the total decay from the 
plateau f A to the equilibrium value zero, i.e., 4>a{t~a) — 
/a/20. Within the range of validity of the second scaling 
law, Eq. (|l9|), the ^4-dependent scales are coupled in the 
sense that ta — TAt' a - Here ta is an ^4-specific control- 
parameter-independent factor determined by 4>a (ta) = 
/^/20. Applying these results to the dipole correlator, 
one gets C\(t) = C\{t) for \a\ -C 1 and t& <?C t where 
t = t/t' a . For the a-scale factors one finds fi = 18.8 (8.73, 
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2.66) for C = 1.0 (0.8, 0.6). The description of the as- 
process for elongation parameter £ = 0.8 is demonstrated 
in Fig. 8 of Ref. ||. 

Using the second scaling law for the tagged-particle- 
density correlator, one gets from Eq. (|5a|) the second 
scaling law for the a-process of the MSD: 

A s (t) = A s (t), i = i/C | a | < 1, t a < t. (29) 

An a-relaxation time ts shall be defined by that time, 
where the diffusion asymptote D$ t is reached within 5%, 
i.e., As(ts) — 1.05 Dsts- One gets Ds — Ds/t' a and 
T s = Tst' a , where Ds and fs are to be determined from 
lim^ A s (t)/t = D s and A s {f s )/fs = 1.05D S . One 
finds Ds = 0.0171 and fs — 11.6. One gets in particular 
Ds cx l/t' a cx (tp c ~ (p) 7 , and Fig. || demonstrates how 
the diffusivity approaches this power-law asymptote for 
ip increasing to ip c . The asymptotic description of A<j(i) 
by Eq. @ is demonstrated in Fig. 9 of Ref. |§§]. 

For times of order t a , the relative corrections to the 
first scaling law are of order y/\ a ]. For times of order 
t'p, the corrections to the second scaling law are of order 
| a |. Therefore, the second scaling law holds for larger 
separations | a | than the first one |22j| . For example, even 

for the large value e = —0.1, D 1 ^ 1 differs from the linear 
asymptote by only 15%, as is demonstrated by the n = 3 
result in Fig. [| The corrections to Eq. ( p9| ) increase if 
t decreases towards t a . But for t « t a , the description 
by the first scaling law becomes valid, which provides 
the leading corrections to Eq. (|29|). The descriptions by 
the two scaling laws overlap. Together, they provide a 
complete description of the dynamics for t > t*. This 
holds provided | e | is small enough, as is demonstrated in 
Fig. | for e = -10~ 3 for A s (t) and A c (t) for £ = 0.8, and 
for A^(t) for 4 values of £. From the analogous figure for 
e = — 10 -2 one concludes that the scaling-law description 
accounts for the MSD quantitatively for t > t* and | e | 
< 0.01. For larger distance parameters, corrections to 
scaling become visible. 

D. Rotation-translation-coupling effects 

It might be adequate to start the discussion of rotation- 
translation-coupling effects by two side remarks. Firstly, 
it was shown that the correlators for the dipole and the 
quadrupole dynamics for £ = 0.8 are in semi-quantitative 
agreement with the experimental data for propylene car- 
bonate |3jJ. Thus, the results to be discussed for £ > 0.6 
can be considered as relevant for the interpretation of 
glass-forming van-der-Waals liquids. Secondly, the sys- 
tem under study exhibits two glass phases for ip > <p c . 
There is a critical elongation £ c so that for £ > £ c cor- 
relations of the molecule's axis arrest for long times as 
do all other correlations of variables characterizing the 
structure. In particular C\{t — > oo) = /i > 0. But 
for £ < £ c , dipole correlations exhibit ergodic behavior, 



i.e. Ci(t — > oo) = 0. Precursor effects of this glass- 
glass transition at £ c disturb the standard transition sce- 
nario p5| , p5| . This is the reason why A^(i) for ( = 0.4 is 
remarkably different from A^(£) for the other £ shown 
in Fig. |^; for £ = 0.4, the elongation is too close to 
C (: = 0.380. For example, for £ = 0.4, the corrections 
to von Schweidler's law, which shift the dashed lines on 
to the full ones for t a < t < t' a , are negative, while they 
are positive for the other cases. 

The glass-transition of the HSS is driven by density 
fluctuations with wave numbers near the first-peak po- 
sition of the structure factor S q close to q — 7.0/d. For 
the scale factor for the a-process of these fluctuations, 
denned by 4> q (f q ) — /q/20, one gets 77. = 6.0. For 
the a-relaxation of the tagged-particle-density correla- 
tions for the same wave vector, one gets a similar num- 
ber fy = 5.0. One concludes that the a-scale factor 
fs = Ts/t' a = 11.6 for the approach of Ag(i) to the dif- 
fusion limit is in the range within which relevant density 
fluctuations decay to zero. The corresponding interval 
for the crossover from the end of the von-Schweidler de- 
cay to the beginning of diffusion is 2.2 decades. This is 
shown by the lowest curve in Fig. [)] and should be consid- 
ered as the normal behavior for the density-fluctuation 
dynamics in simple systems. The MSD for the molecule's 
center behaves quite similarly; Ac(i) for £ = 0.8 is nearly 
indistinguishable from As(t). However, Fig. || demon- 
strates also that A^(t) behaves differently. For ( = 1.0 
(0.8, 0.6), the above specified crossover intervals are 3.3 
(3.1, 2.8) decades wide. This means, the crossover in- 
tervals of Ayi(t) are larger than those of As(t) or Ac(t) 
by factors of about 13 (7, 4), and so are the a-scales ta 
compared to T5 or re- The reason is the reorientational 
contribution to the MSD of the constituent atom, i.e., 
the second term on the r.h.s. of Eq. (||). There are two 
effects. First, for t/t' a ^S> ft, the dipole correlations are 
decayed to zero, i.e., the molecule's axis is distributed 
on the unit sphere with a constant probability density. 
Therefore, there is the positive offset X — (£d) 2 /12 be- 
tween the two functions referring to the atom and the 
center. A^t) — Ac{t) = X for t = t/t'^ > fi as is demon- 
strated in the linear representation of the a-process mas- 
ter functions in Fig. 10. The slow decay of the relative 
offset X/Dct explains that the ratio ta/tc is larger than 
unity and this the more the larger £ is. Second, steric hin- 
drance for reorientations increases with £. For £ > 0.8, 
the a-scale f% for dipole relaxations becomes comparable 
to or larger than fc- Therefore, there appears a interval 
1 < t/t' a < fi, after the end of the von-Schweidler-law 
regime and before the beginning of the diffusion regime, 
where the A J 4(t)-versus-t diagram, in contrast to the 
Ac(i)-versus-£ one, exhibits a curvature. This is shown 
for £ = 0.8 and 1.0 in Fig. [h]. 

If a hard sphere gets expanded to a dumbbell with 
a small elongation, the molecule's center gets restricted 
more tightly in its cage. Therefore, provided £ is small, 
the localization length r c c is smaller than r c s and it de- 
creases with increasing £. For £ 1.0, the localization of 
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one constituent atom of the molecule restricts the motion 
of its partner. Therefore, Tq is smaller than r c s also for 
large elongations. Counter intuitively, the theory does 
not lead to a monotonic interpolation between the speci- 
fied limits. The r£. 2 -versus-£ diagram in Fig. |ll| exhibits 
an oscillation and exceeds r c s by about 10% for £ near 
0.6. 

For ^ < Cc one gets C\{t — > oo) = Q. If one neglects 
the variation of the r£, 2 -versus-£ curve, Eq. leads 
to r\ 2 = r% 2 + (CO 2 / 12 - Th e dotted line in Fig. 
shows that this formula explains the increase of r c A for 
£ increasing up to £ c . For £ increasing above £ C) the 
decrease of 1 — /i and of Tq explains the decrease of r c A . 

It seems plausible that the number of paths for the 
molecule through the system's configuration space de- 
creases with increasing elongation. Therefore, one might 
expect that the diffusivity D A — Dq is smaller than Ds 
and decreases with increasing £. Figure |ll| confirms this 
expectation for small and for large elongations. But, 
strangely, the calculated diffusivity is not a monotone 
function of £; and for £ w 0.6, Dq exceeds Dg by about 
30%. 



IV. CONCLUSIONS 

MCT predicts that the long-time parts of the density- 
correlation functions <f) q (t) of the solvent and <j>q{t) of the 
solute deal with relaxation in the sense that they are su- 
perpositions of Debye-relaxation processes exp[— j(t/to)]. 
They deal with structural dynamics in the sense that 
the weight functions p q {^) and pg(j) for these super- 
positions are uniquely determined by the equilibrium 
structure. The corresponding dynamics is therefore re- 
ferred to as structural relaxation. Within the structural- 
relaxation regime, the subtleties of the normal-liquid dy- 
namics merely enter via the value for an over-all time 
scale to, which is defined with the aid of the critical de- 
cay law, Eq. (|l6|) . The structural-relaxation interval was 
estimated for the hard-s phere system (HSS) under study 
as t > 20t (Sec. [II B , Fig. |j). The specified proper- 
ties of the MCT dynamics are asymptotic ones, valid for 
large times near the transition point. The estimation 
of the range of validity depends on the accuracy level 
required and also on the particular function under dis- 
cussion Differentiation with respect to time, as is consid- 
ered in Eq. (Q) for the velocity correlator Ks(t), enhances 
oscillation features, i.e. deviations from relaxation be- 
havior. Figure ^ demonstrates indeed, that one would 
estimate the structural relaxation interval for Ks(t) as 
t > 40£o, i-e. somewhat more restrictive than for the 
MSD. 

For a normal liquid, the velocity correlator is positive 
for short times. It is also positive for times exceeding 
a cross-over time th, where it exhibits a hydrodynamic 
long-time tail proportional to i -3 / 2 Q]. But, upon ap- 
proaching the glass-transition point, the regime for hy- 



drodynamics shrinks to lower frequencies. Therefore, the 
time th diverges if the packing fraction ip approaches the 
critical value ip c . The cage effect causes Ks(t) to be neg- 
ative for t > 20£o- Even more, the correlation of the 
particle velocity at time t with the opposite of its initial 
value, say = (vs(t) ■ [—Vs(t = 0)]), is proportional to 
a completely monotone function F(t/to), Eq. (p5|). The 
function $ represents a relaxation process in the sense of 
Eq. (23b), where the weight function ^0(7) is given by the 
structural function. These results for the velocity corre- 
lator for t > 20to express concisely the essence of glassy 
dynamics, namely the relaxation caused by the cage ef- 
fect and determined by the Boltzmann factors for the 
equilibrium structure. 

The MSD increases with time according to the charac- 
teristic two-step pattern of the MCT-transition scenario. 
There is a stretched approach towards a plateau value 
rg 2 followed by arrest at the square of the localization 
length r| < r L s 2 within the glass or by the stretched 
start of the a-process within the liquid, Fig. 0. The sen- 
sitive density dependence of the MSD near the plateau is 
described by the universal formulas for the first scaling 
law, as is demonstrated in Fig. ^. The same holds for 
the dipole correlator C\(t). This is shown in Fig. || by 
scaling C\{t) for three values of the molecule's elongation 
on the same functions \J\a\g-{t/t a ). However, there is 
a large time interval at the beginning of the structural- 
relaxation regime, 2Qt < t < t* w 800io, where the uni- 
versal leading-order-asymptotic formulas do not describe 
the MSD, as is shown in Figs. ^ and |7|. The theory pre- 
dicts similar results for all systems with a structure sim- 
ilar to that of the HSS, i.e. for all van-der-Waals liquids. 
Obviously, it would be worthwhile to test by experiment 
or molecular-dynamics simulation whether the specified 
prediction is correct, in particular, whether MCT can re- 
produce properly the structural relaxation of the MSD 
outside the scaling-law regime. 

The beginning of the a-process of the MSD, i.e. the 
increase of A x (t) above the plateau r^ 2 , is described by 
von Schweidler's law, Eq. (|2q). It is exhibited in Fig. ^ 
for t > t a by the dashed lines. The a-process terminates 
in the diffusion law for long times, exhibited in Fig. |9|by 
the straight dashed-dotted lines. The a-process follows 
well the second scaling law, Eq. (|29|), which is presented 
by the dotted lines. The cross-over interval from the end 
of the von-Schweidler-law description (indicated by the 
open diamonds) to the beginning of the diffusion pro- 
cess (indicated by the filled diamonds) for the MSD of 
an atom is about 2 decades wide, as shown by the low- 
est curve in Fig. Q The same is true for the MSD of 
the molecule's center. But, the cross-over interval for the 
MSD of the constituent atom of the molecule is much 
larger due to the rotation-translation coupling. The ex- 
pansion is more than an order of magnitude for £ = 1.0. 
It would be valuable for an assessment of the MCT for 
molecular systems if this prediction could be tested by 
molecular-dynamics simulation. 







Bennemann et al. [[36| have determined by molecular- 
dynamics simulation the MSD for the monomer and for 
the center of a polymer for a model of a glassy polymer 
melt. They identified for their data the a-process in the 
sense of MCT and interpreted it consistently with the 
universal asymptotic formulas. They observed the ex- 
pansion of the crossover interval for the monomer MSD 
relative to that for the center similar to what is shown in 
Fig. U for C = 1-0, and they attribute this expansion to 
the glassy dynamics of the Rouse modes. These modes 
for the polymer's internal degrees of freedom arc thus 
identified as part of the a-process [[36|. Figures ^ and 
|l0| corroborate their conclusions; albeit for the simplest 
polymer model only, namely a rigid hard-sphere dimer. 

Let us consider a one-dimensional model, where two 
hard spheres are restricted to move in a finite interval. 
One calculates easily the localization length r$ of one of 
the spheres, assuming the other sphere to move freely 
between the wall and its partner. Similarly, one can cal- 
culate the localization length rc of the center of a dumb- 
bell built by the two spheres. One finds rc > rs, i.e., 
the freely moving partner provides a stronger hindrance 
for the motion in the cage than the bonded one. Pos- 
sibly, the result r c c > r c s shown in Fig. 11 for £ ~ 0.6 
is the analogue of this phenomenon for one-dimensional 
localization. To corroborate this reasoning, Fig. [l^ ex- 
hibits as full lines the isotropic part (r) of the solute- 
solvent pair distribution function and the corresponding 
structure factor So(q) for ip — (p c for a dumbbell with 
C = 0.6 |^9|. The dashed lines exhibit the results for 
a sphere with diameter d c g = 1.215 d chosen such that 
its volume agrees with that of the dumbbell. The well- 
known excluded-volume effects for this effective spheri- 
cal solute are larger than those for the tagged particle 
of the HSS. In particular, the structure factor peak is 
higher. Therefore, the localization length r^ s = 0.0588 d 
for the effective spherical solute is much shorter than the 
localization length r c s = 0.0746 c? for the tagged particle 
of the HSS. The ratio Tg ff /r| decreases with increasing 
The pair distribution function of the effective sphere 
agrees reasonably well with go (r) for the dumbbell for 
r > 1.5 d. Therefore, the two structure factors are close 
to each other on the wings of the Sq (q)-peaks and also for 
small q. However, packing cannot be done so efficiently 
around the dumbbell than around a sphere of the same 
volume. Therefore, the pair distribution function for the 
latter is much bigger for r near d e g than (r) for the 
dumbbell. This depletion effect leads to a reduction of 
| Sg(q) | at the peak position and also for larger wave 
numbers. Hence, the depletion reduces the magnitude 
of the mode-coupling coefficients. As a result, the local- 
ization length Tq = 0.0863 d, calculated by using MCT 
equations for a simple system with only the isotropic part 
of the dumbbell structure factor, is not only larger than 
r£ ff , but it is even larger than r%. The estimated 7q is 
rather close to the true localization length r c c = 0.0808 d 
for the dumbbell's center. Only if £ becomes even larger, 
mode-couplings to the non-isotropic parts of the density 



fluctuations take over and compensate for the depletion 
effects. This leads to the decrease of r c with increasing 
£ for £ > 0.6, as shown in Fig. [ll]. 

The above discussed depletion effect means, that the 
dumbbell is surrounded by a liquid of lower averaged 
density than specified by the packing fraction. Hence, 
the dumbbell moves in a complex, which effectively has 
a larger distance from the glass-transition point than is 
specified by {ip c — ip)/ip c . Such lubrication phenomenon 
for the dumbbell in the complex might explain the in- 
crease of the diffusivity exhibited in Fig. 11 for £ ~ 0.6. 
Diffusion in a dense liquid is a collective phenomenon. A 
particle can move over a distance comparable to its di- 
ameter only, if one of its neighbors moves. This neighbor 
can move only if its neighbor moves and so on. On the 
average, a flow pattern will be built like one knows it for 
the motion of a sphere in an incompressible ideal liquid. 
It was one of the original motivations for the formula- 
tion of the MCT-equation of motion to approximately 
treat the indicated back-flow phenomenon. Molecular- 
dynamics simulation has given evidence that a typical 
event underlying the mentioned averaged back-flow pat- 
tern might consist of quasi-one-dimensional motions of 
clusters of particles p7| . It seems possible, that the bond- 
ing of two atoms to a molecule stabilizes the cluster, pro- 
vided the elongation £g? is smaller than the diameter of 
the ring formed by the moving cluster. If this was the 
case, the diffusivity for the molecule should be larger than 
that for a tagged atom. Hence, it is not obvious that the 
non-monotone variation of the diffusivity with changes 
of C, which is shown in Fig. [yj is a mere artifact of the 
approximations underlying the presented theory. 
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APPENDIX A: MODE-COUPLING-THEORY 
MODELS FOR COLLOIDS 

In this appendix, it will be shown how the MCT formu- 
las for simple systems can be specialized to models for the 
dynamics of colloids. Within this frame, a representation 
of the velocity correlator in terms of a completely mono- 
tone function F col (t) shall be derived. To begin, let us re- 
member the definition of Laplace transforms of functions 
of time, say F(t), to functions of the complex frequency s: 
F(s) — J Q dt exp(-st) F(t). Applying this transforma- 
tion to 4> q (t), M™ s (t) and m q (t), one can rewrite Eq. (^) 
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together with the initial conditions <p q (t = 0) = 1 and 
dt<t> q {t = 0) = as a double-fraction for 4> q (s): 

= 1 / {s + Q 2 q I [s + M™z(s) + Q 2 q m q (s)}}. (Al) 

The analogous result can be obtained for the tagged- 
particle correlator 

=!/{*+ (^s) 2 / [* + Mf rcg ( S ) + (qv s ) 2 m s q (s)}}. 

(A2) 

Equation (Q) can be rewritten as a relation between the 
Laplace transforms Ag(s) andifs(s) of Ag(t) and Ks(t), 
respectively: 



K s (s) = 3s 2 A s (s). 



(A3) 



The small-g-expansion of Eq. (A2) leads to 



K s (s) = 3v| / [a -I- M s < reg ( S ) + „| ms(s)}. (A4) 



Here M S ' rog (s) = lim 
+0 q 2 m 



q ^ M^(s), 



and ms{s) 



lim„^n q 2 m^i (s) is the transform of ms(t) from 



Eq. @) @. 

In a colloid, there is a contribution to the fluctuat- 
ing force due to the interaction of the particles with 
the suspending liquid. This leads to contributions to 
the rates M q ° s (s) and M q ,rcs (s) which are large com- 
pared to the scale for the frequency s one wants to con- 
sider. Therefore, it is the conventional first approxi- 
mation step of a theory for colloid dynamics to coarse 
grain correlation functions on time scales of the order 
of 1/Af 9 r °s(s) and 1/M^' reg (s). For the correlators con- 



M g rc s(s) in Eq. (JAJ]), 



sidered here, this means that s 
s + Mf' reg (s) in Eq. (|Aj), and s 

are replaced by v q = M r q c s(s = 0), if = M^ TC ^(s = 0), 



M S ' reg (s) in Eq. (IAJ) 



; s - 



(s = 0) = lim 9 ^o v q ! respectively. Let 
us indicate the various functions, obtained by this spe- 
cialization of MCT, by superscripts "col" . The equations 
of motion, obta ine d by t he b ack transformation of the so 
modified Eqs. (O) and @, are 



mf(t)=^[C 1 (*)], mf» col (t)=^| 



,S,col 



(*),<(*)]• 

(A7) 



Equations (|A5|)-( A7), together with the initial conditions 
cj) q °\t = 0) = 1 and </>f col (i = 0) = 1 define a unique 
solution with all the general properties of correlation 
functions |52j. Since the coarse graining has altered the 
s — > oo asymptote of the corr elator tra nsforms relative 
to that exhibited by Eqs. ( |A1[ ) and ( |A2[ ), the short-time 
behavior is altered, namely 4> q ol (t — * 0) = 1 — \t\/r q + ■■■■, 

c/)^ col {t — > 0) = 1 — \t\/r q -\ . Since s was dropped 

in the denominator of Eq. (A4), the velocity-correlator 



Laplace transform for the colloid dynamics does not tend 
to zero for large s; rather Kg ol (s — > oo) = 3Dq , with the 
abbreviation Dq = v 2 s /vq. This suggests to rewrite the 
coarse-grained velocity correlator as 



D S2pco\ 



F 



( S )/3 

m^( S ) /[1 + D. 



(s), 



col 



00], 



-.col 



(t)=^M' rol (*), <(*)]• 



(A8a) 
(A8b) 
(A8c) 



The solutions of Eqs. (|A5|)-( A7) are completely mono- 
tone @. Kernel m c s °\t) is a combination with positive 
coefficients of products of completely monotone functions 
because of Eqs. ( A8c ) and (H) , and therefore it is com- 
pletely monotone as well J33f . Because of Bernstein's 
theorem, one can write (s) = J Q [s + 7] _ V(7)<^7 
with a non- negative distribution p(j)- Therefore, there 
hold the following four properties |3j| for m™'(s): (i) it 
is holomorphic for all complex s, except for negative real 
numbers, (ii) m c g l (s)* = m c c? l (s*), (iii) m c s ? l (s) — > for 
Res — > oo, and (iv) Imro™ (s) < for Ims > 0. One 
concludes furthermore that 1 + Dq m c g l (s) ^ for all s 



except, possibly, negative real values. From Eq. (A8b) 
one concludes that also F col (s) exhibits the properties 
(i)-(iv). Hence, F col (s) is the Laplace transform of a 
completely monotone function F co1 (t) (Rcf . |33| , ch apte r 
5, theorem 2.6). The back transformation of Eq 
leads to the desired representation 



(A8a 



K c s °\t)/3 = D^5(t)-DpF co \t), 



(A9) 



T q d t <tf q \t)+Cjf q \t) 
ft 



+ / dt'm q \t-t')d t ,<t> q \t') = Q, (A5) 
Jo 

ft 

dt' mf' co \t - t') d t ^ s q co \t') = 0. (A6) 

o 

Here T q = v q j^l? q and r q — v q /{qvs) 2 are times charac- 
terizing conventional colloid dynamics. The interaction 
potentials between the particles, which cause the cage ef- 
fect, are not altered by the introduction of the solvent. 
Therefore, the expressions for the mode-coupling contri- 
butions to the kernels, Eqs. ^ and (Q), keep their func- 
tional form 



where F co \t) obeys Eqs. (§|). 

Let us add, that the back transformation of Eq. ( A8b| ) 
leads to an equation of motion determining F col (t) from 
'(*): 



the kernel m™ 1 



F co \t) = m c s ol (t) - D$ / dt 1 m c s ?\t - t') F co \t'). 
Jo 

(A10) 

Integratin g Eq. ( |A9| ) twice with respect to the time and 
using Eq. (A10), one gets an equation of motion for the 
MSD of the colloid 



t- / dt'm c s oi (t-t') A c s ol (t') 
o 



(All) 
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This result was derived originally along a different 
route |§]. It implies lim. t _> A^ ol (t)/t = D$. 
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FIG. 1. Normal- 

ized velocity-correlation functions Ks(t)/(3vg) for a tagged 
particle of the hard-sphere system (HSS). The dotted lines 
with label c refer to the critical packing fraction ip c w 0.516. 
The full lines with labels n = 2 and 3 are calculated for dis- 
tance parameters e = (tp — ip c )/ip c — ^fl0~ n ^ 3 for the liquid 
(e < 0) and for the glass (e > 0), respectively. Here and in 
some of the following figures an arrow marks the time 20to, 
where to = ic? = 0.0236(d/t;) is the time scale from Eqs. ( p^ ) 
and ( |2lj ) for the critical decay. The units of length and of 
time are chosen here and in the following figures so that the 
particle diameter d and the thermal velocities v = vs are 
unity. 
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FIG. 2. Double logarithmic presentation of 
A s (t) = Sr 2 s (t)/6 for the mean-squared displacement (MSD) 
8r%{t) for a tagged particle of the HSS. The dotted line with 
label c refers to the critical packing fraction tp c and the full 
ones to e = ±10~ n//3 . The straight dashed line with slope 2 
exhibits the ballistic asymptote (vst) 2 /2. The straight dot- 
ted lines with slope 1 exhibit the long-time asymptotes Ds t 
of the two liquid curves for n = 1 and 9. The horizontal 
line marks the square of the localization length at ip — ip c : 
r c s 2 = 0.00557 d 2 . 
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FIG. 3. The diamonds are the values for D]^ 1 with the 
HSS exponent 7 = 2.46 for the tagged particle diffusivities 
Ds determined for the liquid curves in Fig. ^. The straight 
line is the function T(tp c — (p), ip < (p c , with V chosen so that 
the line goes through the data point for n = 9. The crosses 
exhibit the square of the localization length, r%, determined 



for the glass curves in Fig. 
asymptotic law r c s 2 — hs \ 



2j. The full line exhibits the leading 
f/VT=\, Eq. (§§). 
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FIG. 4. Double logarithmic presentation of the MSD for a 
tagged particle of the HSS as function of the reduced time t/to 
for distance parameters e — ±10~ x ; x — 1,2, 3. The full lines 
reproduce the results from Fig. ^ with labels n = 3, 6, 9 and 
to = tf from Eq. @. The dashed lines are the corresponding 
results for the colloid model specified in the text J23] with 
to ~ to from Eq. (^) . The curves for x = 2 and 3 are shifted 
by two and four decades, respectively, to the right to avoid 
overcrowding. The full dots and squares ma rk the time scales 
tcr and t' a , respectively, defined in Eqs. (17c) and (18b). 
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FIG. 5. Rescaled velocity corre- 

lators -F(t) = toKs(t)/3 = tod?A s {t) as function of the 
rescaled time t/to for a tagged particle of the HSS at the crit- 
ical packing fraction ip — tp c . The full line refers to the model 
for a Hamiltonian dynamics with to = to from Eq. and 
the dashed one to the colloid model defined in the text with 
t =1% fromEq. @. The full line is a rescaling of the dotted 
curves in Fig. |l| where values for t < 0.2, i.e. t/to < 8.5, are 
not reproduced. 
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FIG. 6. A s (t) for a tagged particle of the HSS 
at the transition point ip — tp c (full line), the 
leading-asymptotic expansion, Eq. ( prj ) (dashed line), 
and the leading-plus-next-to-leading-asymptotic expansion 
As(t) = r c s 2 - hs(to/t) a + k s (to/t) 2a (dashed-dotted line). 
The diamond and circle mark the times t* = 18.7 = 792to 
and t** — 1.55 = 65.7to where the full line differs by 5% from 
the dashed and dashed-dotted line, respectively. The horizon- 
tal line marks the long-time asymptote r| 2 . The dotted line 
exhibits the ballistic asymptote ^(vst) 2 . 
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FIG. 7. The full lines show As(t) for packing fractions ip 
given by e = (tp — <pc)/<fic = —W~ x , x = 2, 3, 4. The dashed 
lines are the first-scaling-law descriptions by Eq. (|2t[) . The 
diamonds mark the points where the dashed lines differ from 
the full ones by 5%. Within these intervals, As(t) varies be- 
tween about 0.0043 and 0.0086 as is indicated by the shaded 
bar. The crosses exhibit the von Schweidler law, Eq. 
The filled circles and squa res m ark t he t imes t a and t' a , re- 
spectively, defined in Eqs. (17c) and (18b). 
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FIG. 8. The full lines are dipole correlators Ci(t) rescaled 
to Ci(t) — [Ci(t) — fi]/hi for the elongation parameters ( 
= 0.6, 0.8 and 1.0 (from top to bottom) where /f = 0.769, 
0.905, 0.955 and hi = 0.46, 0.19, 0.09, respectively. The 
distance parameters are e = — 10~~ x , x — 2, 3, 4 (compare 
text). The filled circles and squares mark the times t a and 
t' a , respectively, for the three packing fractions. The dashed 
lines exhibit the first-scaling-law asymptotes yj\ a \g~(t/ta). 
The dotted line is the MSD for a tagged particle for x — 3 
rescaled to A s (t) = [r c s 2 - A s (t)]/h s . 



16 




FIG. 9. As(t) and A c (t) for < = 0.8, and A A (t) for C 
= 0.4, 0.6, 0.8, 1.0 (full lines, from bottom to top). Suc- 
cessive curves are shifted upwards by two decades to avoid 
overcrowding. The distance parameter is e = — 10~ 3 and 
the corresponding times t a and t'^ are marked by filled cir- 
cles and squares, respectively. The dashed lines are the 
first-scaling-law asymptotes, Eq. (p7[). The open diamonds 
mark the points where the dashed lines differ from the full 
ones by 5%. The straight dashed-dotted lines exhibit the 
diffusion asymptotes Dst, Dct and and the filled dia- 

monds mark the position where these differ from the full lines 
by 5%. The dotted lines, which coincide with the full ones for 
t > 10 4 , exhibit the second-scaling-law asymptotes, Eq. 
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FIG. 10. a-relaxation master functions A as function of 
the rescaled time t = tjt' a for ( — 1.0, 0.8 and 0.6 (from top to 
bottom). The curves for £ = 0.8 (1.0) are shifted upward by 
0.4 (0.8) in order to avoid overcrowding. Curves with labels 
C and A refer to the MSD for the molecule's center and for 
the constituent atom, respectively. The filled diamonds mark 
the time fc, where the diffusion asymptote is reached within 
5%: Ac(tc) = 1.05-Dcfc- The corresponding time ta for the 
constituent atom for £ = 1.0 (0.8, 0.6) is f A = 140 (72, 39). 
The triangles mark the time fi, where the dipole correlator 
has completed 95% of its a-decay: Ca(ri)//i = 0.05. 
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FIG. 11. Squares and triangles are the localization lengths 
squared at the transition point if = tp c as function of the 
molecule's elongation £ for the constituent atom (r^ 2 ) and 
the center (r^ 2 ), respectively. The circles are the diffusivity 
Da ~ Dc of the molecule calculated for e = —0.01. The lines 
connecting the symbols are guides to the eye. The horizontal 
dashed line is a common one indicating both the values of the 
square of the localization length r c s 2 (on the right scale) and 
of the diffusivity Ds (on the left scale) for a tagged particle 
of the HSS. The dotted line is the function r c s 2 + (<d) 2 / '12 
discussed in the text. The inset exhibits 1 — /i as function of 
£, where /i = Ci(£ — * oo) is the long-time limit of the dipole 
correlator Ci(t). The arrow indicates the critical elongation 
Cc = 0.380 for a glass-glass transition. 



to 




1.0 



2.0 
r 



3.0 




FIG. 12. Pair distribution g (r) as function of the dis- 
tance r between solvent particles and the center of the so- 
lute (upper panel) and the corresponding solute-solvent struc- 
ture factor S s (q) (lower panel) calculated within the Per- 
cus-Yevick theory for the critical packing fraction ip c = 0.516. 
The dotted lines refer to a tagged sphere of the same diame- 
ter, ds = d, as the one of the solvent particles. The dashed 
lines refer to a spherical solute of diameter d e g = 1.215 d (see 
the text). The full lines exhibit the isotropic part of the dis- 
tribution g"o(r) and the structure factor 5*o^(g) for a dumbbell 
consisting of two fused hard spheres of diameter d with the 
elongation £ = 0.6 p9|. 
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